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Rigorous Electromagnetic Analysis of Uncooled Microbolometer 


I. Problem Statement 

In the last decade, thermal infrared detectors that are operated near room tempaature have made large advances in performance, 
so that they are now used in, or being considered for, a wide range of commercial and military applications. Thermal detectors 
work by converting incidait infrared radiatirm to heat, thereby inaeasing the temperature of the detector. The change in 
temperature is typically measured as a change in detector resistance (for bolometers) or in the spontaneous polarization across 
the detector (for ferroelectrics). 

For best performance, the detector should have the highest possible absorption, lowest possible thermal conductance between the 
detector and its surroundings, and a thermal time constant roughly equal to the frame rate. These requiremaits are sometimes 
incompatible. For example, if the thermal conductance is reduced without also reducing the detector mass, the thamal time 
constant becomes too long. Therefore, fiirther improvements in sensitivity will require reductions in thermal conductance and 
mass while maintaining high abs(»ption. 

There is also a trend to reduce pixel size in order to improve image resoluticm. However, as pixel sizes shrink to only a couple of 
optical wavelengths of the incident radiation, it is unclear how diffraction will effect optical absorption. To investigate small 
pixel effects it is necessary to employ a rigorous three-dimensional electromagnetic model. In this iffoject, we used the finite- 
difference time-domain (TOTD) mettod to compute absorptioi in finite size pixels placed within an infinitely periodic imaging 
array. The model was used to determine the fractional absorption as a function of wavelength, determine the effect of fill-factor 
and pixel size on absorpticm, and was used to compute an image performance model function. 

II. Physical EteTECTOR Model AND Parameters 

Figure 1 shows a typical thermal detector. The design uses a quarter-wavelength cavity [1,2] to achieve high absorption in a 
thin, low-mass structure. The ideal form of the quarter-wavelength structure is an infinite conducting film with sheet resistivity 
equal to the free-space impedance Z© = llOit £2, spaced one-fourth of a wavelength (k/4) above a perfectly reflecting surface. 
Such a structure has a theoretical absorption of 100 percent at the design wavelength. 


Inddeit Radiaticxi 



Support/Isdation 
Leg 


Figure 1: Typical thermal detector structure 
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In an actual uncooled focal plane array (FPA), the metallized surface of the readout integrated circuit (ROIC) acts as the 
reflecting surface, and the bolometric or ferroelectric material is formed into a thin sheet spaced about 3 pm above the ROIC to 
form the quarter-wavelength cavity. Typical pixel pitch in current uncooled focal plane arrays (FPAs) is 50 pm, with a decrease 
to 25 pm anticipated in the next generation. Since this is only two to three wavelengths of the 8-14 pm light in the long 
wavelength infrared (LWIR) band, significant small-pixel effects are expected. 

The supporting legs can also affect optical absorption in small pixels. This is especially true when the legs are positioned under 
the absorbing elemait in an attempt to inaease the fill factor. In that case, the legs are inside the optical cavity. 

A one-dimensional transfer matrix method was previously used to compute the absorption properties of large (infinite) pixels [3]. 
That model was used to predict the variation of absorption with cavity thickness and absorber conductance. Use of FDTD also 
permits calculation of small-pixel and leg effects. 

In this study, we were primarily interested in modeling the generic effects of pixel pitch and fill factor, so we have approximated 
the absorbing detector layer as a two-dimensional conducting sheet spaced one quarter wavelength over a perfectly conducting 
backplane. This suppressed effects related to the finite thickness of the detector as well as the supporting legs. The previous one¬ 
dimensional modeling [3] showed that this is a reasonable approximation, although the detailed layer structure of the absorbing 
layer could easily be included in FDTD calculations for specific structures. 

III. FDTD Model 

The FDTD method was first introduced by Yee in his 1966 seminal paper [4]. The method was later refined by Taflove and 
others [5,6,7] and is now one of the most widely used numerical techniques for solving electromagnetic (EM) problems. The 
method is extremely general in the materials and geometries it can analyze. Structures that contain inhmiogenous, lossy or even 
anisotropic material properties can be easily handled 

The FDTD method derives its name fi'om a direct finite difference approximation to Maxwell’s time-dependent curl equations. 

-^ = --Vx£,-^ = -Vxff--£ (1) 

dt 11 dt e e 

Yee’s basic scheme is to use a series of finite-difference approximations to transform Maxwell’s equations into a system of six 
algebraic equations. These equatims can then be solved using a simple “leap-frog” algorithm. For a detailed description and 
derivation of the full 3D FDTD equations the reader is referred to [5]. 

The unknown field distribution over a finite space is calculated by applying the six FDTD equaticns to a volumetrically sampled 
grid of cells called Yee cells. Each edge of a Yee cell may be assigned independent electric^ properties, which allows one to 
model complex objects that consist of lossy, inhomogeneous or anisotropic materials. Figure 2 illustrates how the field 
components are assigned to a typical unit Yee cell. 
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Figure 2: Three dimeasional Yee Cell 

After mcorporating a time-dependent incident field, a time-marching algorithm is employed to calculate the unknown field 
distribution. At t=0, a plane-wave source of frequency f is assumed to be turned on. Solving the finite-difference equations at all 
lattice points simulates the propagation of waves from this source. Time stepping continues until the field components at each 
lattice point reach sinusoid^ steady state. Alternatively, non-sinusoidal incident fields that contain an entire band of firequencies 
(e.g., Gaussian-modulated sinusoids) can be applied. When analyzing structures composed of only linear materials the various 
firequencies of the incident field propagate independently through the FDTD solution space. Therefore, application of a Fourier 
Transform to the time-dependent fields at each point in the solution space yields the steady-state field solution throughout the 
entire frequency spectrum of the incident field, hi other words, we can obtain the entire l<mg wavelength response (8-14 pm) of 
a bolometer array in a single FDTD run. 

All calculaticais were performed using custom FDTD code that contains features that are not commonly found in other EM 
codes. These include Ae ability to exploit inherent symmetries found in some of the more common IR detectors. For instance, 
the bolometer model used here possesses a four fold symmetry that when exploited can significantly reduce the computational 
complexity of the problem. They also include the ability to simulate infinite periodic arrays by selecting appropriate boundary 
conditicns (BCs). 


IV. Summary of Most Important Results 
A. Optical Absorption as a Function of Wavelength 

We used the FDTD model to calculate the fractional absorption (total absorbed power divided by total incident power) of infinite 
and finite size pixels throughout the long wave IR band (8-14 pm) 

1) Infinite Size Pixel 

We first compared the FDTD result to previous calculations obtained using a transfer-matrix approach used for infinite pixels 
[3]. The purpose of these calculations was to validate the FDTD model using an independent analytical model and to refine the 
FDTD parameters (cell size, time step, material properties) to be used in the finite-pixel size calculations. The structure 
consisted of a conducting sheet with sheet conductivity of 1/^ , where = 1207tQ is the impedance of free space, spaced 

2.5 pm above a perfectly conducting surface (X/4 cavity at a 10 pm. freespace wavelength). Periodic BCs were used in the 
FDTD runs to simulate an infinite pixel. Such a structure has a theoretical absorption of 100 percent at a wavelength of 10 pm. 
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Figure 3 shows the computed spectral absorption for normally incident light. Note that the FDTD result for the two smaller grid 
sizes is very close to the analytical (transfer matrix) result. Based on this result, we used a grid size of 0.25 pm for the remainder 
of the calculaticHis reported here. 



Wavelength (gm) 


Figure 3: Spectral absorption computer using analytical (transfer matrix) method and using FDTD with three cell 
sizes. 


2) Finite Size Pixel 

A series of pixels were run that consisted of a square detector centered within a larger square unit cell. Periodic BCs were 
applied at the sides so that the detector was implicitly embedded in a periodic array of other identical detectors as in a focal plane 
array (FPA). The fill factor, F, of each pattern was defined as the fraction of area in the unit cell occupied by the detector. Figure 
4 shows the calculated absorption for several finite geometries compared with the infinite-pixel result. We expect the fractional 
absorption, curacy to be approximately linearly dependent on F, i.e., Ofrac = fiF, where tj is a normalized absorption. However, a 
plot of the effective rj versus wavelength. Figure 3, shows that rj is larger in the finite, 25-pm period pixels than in the infinite 
pixel. This implies that absorption decreases more slowly than detector area, i.e. that the effective optical area of the detector is 
somewhat larger than its physical area. 

B. Optical Coupling Efficiency 

To better quantify the expected performance of an uncooled bolometer array we define a quantity, we call the optical coupling 
efficiency. This term integrates the firacticaial absorption over the long wave IR band (8-14 pm) normalized by the black body 
spectrum. Mathematically ^ is written as 


X=iA fan 




X = \Atan 

J BB{X)dX 

X=1 fan 


where BB(X) denotes the black body spectrum. 


( 2 ) 
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Figure 4: Absorption for finite and infinite pixels 


1) Variation with Pixel Size and Fill Factor 

We calculated the optical coupling efficiency of pixels ranging in size from 10-50 pm and area fill factors varying from 0-100%. 
The results are summarized in Figure 5. Without considering diffiacticm effects the coupling efficiency would vary linearly with 
area fill-factor. This is purely a consequence of the linear increase in surface area. However, it is clear from Figure 5, that die 
coupling efficiency variation is far from linear. This is especially true for the smaller 10-20pm pixels. Diffraction effects 
actually improve the coupling efficiency by bending light incidoit outside the physical detector surface into the cavity region 
allowing it to be absorbed after reflection. The optical area of the detector is thus larger than the physical area. For smaller 



Area Fill Factor 

Figure 5. Coupling Efficiency vs Fill Factor for pixel pitches varying from 10-50 pm 
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pixels this inarease in optical area plays a more significant role in the total coupling efficiency (diffraction + geometry) than for 
larger pixels. To obtain a better understanding of this effect we calculated the absorption pattern of pixels ranging in size from 
10-50 pm at a single wavelength (10 pm). Image plots of these absorption patterns are shown in Figure 6. Without diffraction 
effects (Mie would expect a uniform absorption of energy throughout the aitire detector area. For the case of the largest pixel (50 
pm) in Rgure 6 this is nearly the situation except for some increased absorption at the edges and comers of the detector. 
Howevor, as the detector size decreases the absorption pattern becomes hi^ly non-imiform indicating a more significant 
contribution frcm diffractitm. 


10 micron pixel with a 0.5625 %areafiil factor 



Microns 


20 micron pixei with a 0.5625 %area fill factor 



-10 -5 0 5 10 

Microns 



Figure 6: Absorption patterns for incident radiation at a 10p,m wavelength 
C. Image Performance Model Function 


The coupling efficiency shown in Figure 6 includes the geometric effects of a reduced fill factor. The geometric effects are 
included in most imaging system performance models, so we normalize these results to include cmly the coupling efficiency due 
to diffraction. The results are shown in Figure 7. The data from the rigorous models are shown as points. Note that when the 
pixel pitch is reduced to 10 pm with a small fill factor, diffractioi can increase the flux gathering capabilities of the detector by a 
factor of 3 above the geometric estimate of flux gathering capabilities. Even with a reasonable fill factor of 42 percoit, the 10 
pm detector pitch obtains 1.8 times the geometric approximation of flux absorptioi. 

The data points were fit using a three dimensional curve-fitting program, where the curve fit was specified to be simple 
polynomials. A curve fit was obtained with a 0.9844 correlaticm coefficient and is shown as the mesh surface in figure 7. The 
equation for the curve fit approximation of diffraction coupling efficiency was 

=0.957-1.50a-H26.10/^-H4.36a' -39.50/<J" -47.5Sa/S-2.SS3a^ +216.3/-4A14a/ -i-23.94a '/S 

(3) 

where a was the fill factor and 8 was the pixel pitch (micrometers). This equation allows the straightforward calculaticHi of 
diffraction efficiency fm: the longwave miaobolometer as a function of detector pitch and fill factor. This equation can be used 
without the requirement for a rigorous electromagnetic model as long as it is used within the 8 to 14 pm bandwidth and the 
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geometry described above. This efficiency model can be included in infrared imaging system performance estimates to account 
for the effects of diffractim coupling. 



Figure 7: Coupling efficiency due only to diffiactimi (excludes geometric effects) as a function of detector pitch and fill factor. 


To summarize the results, in this (ffoject we modeled the optical absorption of uncooled bolometers using a rigorous three 
dimensional electromagnetic model (FDTD). The model incorporated diffiractitm as well as geometrical effects and quantified 
their relative contribution as a functirai of the pixel’s pitch and area fill factor. The results indicated that diffiaction taids to 
improve the flux gathering capability of a pixel beyond what was predicted using simple geometrical models curraitly employed 
in most imaging system performance models. This effect is more significant as the pixel size is reduced to only a few optical 
wavelengths. We also generated a image performance model function that can be used to include diffiactimi effects in systems 
level modeling. 
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